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Abstract 

We use QCD sum rules for the three point function of a pseudoscalar and 
two nucleonic currents in order to estimate the charge dependence of the pion 
nucleon coupling constant gNNn coming from isospin violation in the strong 
interaction. The effect can be attributed primarily to the difference of the 
quark condensates < uu > and < dd >. For the splitting (g PP n — 9nnno) / 9NNn 
we obtain an interval of 1.2*10 -2 to 3.7*10 -2 , the uncertainties coming mainly 
from the input parameters. The charged pion nucleon coupling is found to be 
the average of g ppno and g nn7Vo . Electromagnetic effects are not included. 
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The effect of isospin violating meson nucleon couplings has recently seen a strong revival 
of interest in the investigation of charge symmetry breaking (CSB) phenomena |1]-§J (for 
a comprehensive review see and references therein). On a microscopical level, isospin 
symmetry is broken by the electromagnetic interaction as well as the mass difference of up 
and down quarks m u ^ m^. It is the aim of this paper to examine the difference between 
the pion nucleon coupling constants g ppn o, g nmz o and g pn7T + using the QCD sum rule method, 
which has been established as a powerful and fruitful technique for describing hadronic 
phenomena at intermediate energies @-^] . Here we will only look at effects which arise from 
isospin breaking in the strong interaction. In the QCD sum rule method this is reflected by 
m u ^ ma as well as by the isospin breaking of the vacuum condensates. Electromagnetic 
effects are not examined. Our work follows the approach of refs. |]7|,|9|JT0| and extends their 
analysis to the isospin violating case. 

We start from the three point function of two nucleonic (Ioffe) |lTJ and one pseudoscalar 
interpolating currents with the appropriate isospin quantum numbers p|,P JT0|jr2]JT3|] , e.g.: 

A NN APuP2,q) = J d 4 Xl d 4 x 2 e ip ^e- ip2X2 (o|T W (2i)/f =1 (0^(^)|o) , (1) 

where i stands for + or and N for proton or neutron, respectively. The expressions for 
the pseudoscalar isovector currents read 

P^ 1 (x) = u(x)i^^u(x) — d(x)i^^d(x) , (2a) 
Plf{x) = V2u{x)i~f 5 d{x), (2b) 

and those for the Ioffe currents are 

r} p (x) = e abc [(w a (x)C7 M u 6 (x)) ~f 5 ^d c {x)] , (3a) 
Vn(x) = e abc [(d^x^^ix)) 7s7^ c (^)] • (3b) 

The momenta pi and p 2 are those of the nucleon, and q = pi — p 2 that of the pion; C = ^7270 
is the charge conjugation matrix. In the following we will only keep terms up to first order 
in isospin violation, i.e. — m u . 



The phenomenological side of the QCD sum rules for the three point functions A are 
obtained by saturating the general expressions for the A's ([1]) with the corresponding nucleon 
and pion intermediate states. In order to connect to hadronic observables we have to know 
the overlap between the pion states and the interpolating fields. The axial Ward identity 

d»A„ a = iq l5 {M,^}q (4) 

gives 

mo(0\uij 5 u - di^dln ) = m w o 2 f w o + 0((m u - m d ) 2 ) (5a) 

V2m (0|m7 5 (i|7r + ) = m 7T + 2 f 7T + + 0((m u - m d ) 2 ) (5b) 

with q = (*), M = m U + ?^f^T 3 and m = ^±IM_ 

Hereby we have used that (m u — md)(0\ui~f5U + <ii75<i|7r ) = 0((m u — md) 2 )- Furthermore 
we can set m n o 2 = m 7T + 2 = m n 2 as well as f n o = f n + = because the differences between 
the charged and the neutral quantities are also of 0((m u — nid) 2 ) [ I3fl . 
We also need the current algebra relation 

m (uu + dd) = (-)m 7r 2 /7r 2 + 0((m u - m d ) 2 ) (6) 

which follows from eq.([5]) and the PCAC relation: 

= m 2 f^\ (7) 

The pion nucleon couplings are defined through the interactions: 

£ppk° = 9ppn0pij5^°P, (8a) 
C nnn o = (-)g n mronij5^°n, (8b) 
C pn7r + = V2g pnK +nijB'T + P< (8c) 

It should be remarked that in our notation all three couplings are positive and have the 
same value in the isospin conserving limit. 
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We then obtain the following expressions for the phenomenological sides of the three 
point functions, eqs.fll]): 

4*"°=*^ —^T^^^^W Mpl5 4 + (9a) 

A - 0=zA " ^^w^W^W Mn75 ' + (9b) 



m^ 2 /* v / 2fi'p„ 7r + 1 1 M p + M n , 
A pnn+ = i X p X n ^-f — — 2 — — ^ 75 f + • • ■ , 9c 



where the A at 's are the overlaps between the Ioffe currents (0) and the corresponding single 
nucleon states. The . . . denote contributions from higher resonance intermediate states and 
the continuum. We will come back to these contributions later. 

By saturating the three point function eq.(|l]) for the neutral current with pseudoscalar 
isovector intermediate states and deriving eqs. (pa|) and (|9b"D we have assumed so far that 
the 7T° mass eigenstate is a pure isovector state. However due to n — rj mixing the correlator 
in eq.(|T]) with the current P^ 1 will pick up a contribution from the \rj > state as wellQ. In 
order to avoid this we have to use a correlator where the pseudoscalar meson current has 
only overlap with the physical |vr >, i.e. the mass eigenstate and not with the \r] >. As 
it has been shown in ref. |l4j] this is possible in lowest order chiral perturbation theory by 
using the linear combination of the SU(3) flavor octet pseudoscalar currents 



P a =3 + 0P a=8 (10) 



where 



T=l 



P a =3 = u(x)i>y 5 u(x) - d(x)i-f 5 d(x) = P i=0 
P _ 1 



u(x)i'j5u(x) + d{x)i^i^d{x) — \f2s(x)i^§s(x) (11) 



rather than the pure isovector current in the correlator eq. ([[]). The 9 denotes the n — 7] 
mixing angle which defines the mass eigenstates \ir > and \rj > in terms of the flavor octet 
eigenstates \n a =3 > and |7r a=8 >: 



1 We are grateful to K. Maltman for pointing this out to us. 
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|vr > = |vr a=3 > +6 |7r a=8 > (12a) 

\V> = K=8 > -0 |7T a=3 > (12b) 

It should be noted in this context that there exists actually a whole family of possible 
choices for interpolating currents involving linear combinations of P a =s and the flavor singlet 
current P a =o, which have no overlap with the r\ but only with the 7r. Our choice fllTf) is the 
appropriate one if one ignores possible mixing to the SU(3) flavor singlet state, i.e. the rf, 
because in this case the current (|ToD is the only choice which has no overlap with the flavor 
singlet state either. 

Furthermore it should be noted that we have neglected all higher pseudoscalar, isovector 
resonances n' , ir" , . . .. In other words we have assumed that pion pole dominance works 
at spacelike q 2 ~ — lGeV 2 , where the three point function method can be applied |7]||. We 
will discuss this point later as well. 

The next step is to perform the operator product expansion (OPE) for the three point 
functions under consideration. Typical diagrams are shown in Fig.[T| Following refs. |7|,|9|JT0 



we keep only terms which are proportional to ^75 and have a \ pole. We identify the residua 
of this pole with one on the phenomenological side, assuming hereby that \q 2 \ ^> m^ 2 , so 
that the pion mass can be neglected in eqs.(|9]). Finally we take p\ 2 = P2 2 = —P 2 in the 
equation of the pole residua and perform a Borel transformation with respect to P 2 . It 
should be noted that the OPE side contains, of course, also terms which do not have a \ 
pole. They will give rise to a form factor, i.e. a q 2 dependence of the pion nucleon couplings 
||13|| , which we do not consider in the present context. 

In our case one can easily convince oneself that up to and including order 4, only the 
diagrams in Figs.|l|(b) and |l|(c), which contain the quark condensates < uu > and < dd > 
contribute. Diagrams containing the gluon condensate < G 2 > (Fig.|I](f)) come in at order 
6, because from dimensional arguments they are proportional to the current quark masses 
m u or rrid, respectively. The mixed condensates < uG ■ ou > and < dG ■ ad > (Fig.|I](g)) 
are genuinely of two orders higher than the quark condensates. Four quark condensates 
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(Fig.|](h)) enter already at order 8. Because reliable values for the isospin breaking of the 
mixed condensates and the four quark condensates are missing, we prefer to stop the OPE 
at order 4 and do not take the higher order power corrections into account. 

Applying the prescription described above one can easily derive the Borel sum rules for 
the three point functions of eq. ([I]): 



7T" 



5 1 jj ' 

- < uu > H — < dd > 

6 6 



+ 







^■7T fir 



5 - 1 Xv 

- < uu > — < dd > 

6 6 



/ 1 r M P - 
V — !L -^ M p (+)g ppn o — e m2 



7T" 



5 J; 1 

- < dd > H — < uu > 

6 6 



m VMV 

r 5 t, 1 

H = - < dd > — < uu > 

L6 6 



rJLv 



(13a) 



A n 2 ^ M n 
m 



)9nnir° [jp ) e ^ 



(13b) 



and 



H-5 



7T" 



- < uu > +- < dd > 
2 2 



M p + M ra / 1 y e mT 



e 



(14) 



It should be noted hereby that the strange quark in the current P a =s (eq.fllOl)) does not 
contribute in the OPE up to that order which we are taking into account. 

Already at this point we see by taking the difference between eq. (|13a|) and eq. ( |13b|) and 



comparing with eq.(14) that up to first order in isospin breaking the charged pion nucleon 
coupling is exactly the arithmetic average of the two neutral pion nucleon couplings, i.e. we 
have: 



9pnir+ 2 \-9ppir 9nnir°\- 



(15) 



which is a simple consequence of the u and d quark contents of the three point functions 
and valid within the approximations considered. 

In order to obtain the splitting between g ppn o and g nnn o we take the sum between eq.( |13a| ) 
and eq.( p.3b| ) and divide by either one of them. Expanding again up to first order in isospin 
breaking, we obtain: 



Sg 



5X r 2 * 



5M 



N 



g.NN J V X N 2 J ' \ M N J ' " \M N J V M 2 J 
Here we have used the following notations for the isospin splittings: 



+ 2 



'5M 



N 



'Ma 2 ' 



'TV 



2 4 e 

~3 7 + 375" 



(16) 



5M N = M n — M,, <5g = # nn7r o — # pp7r o, 5Aat 2 = A n 2 — A p 2 , 

and the average values 

11 1 

M N = -{M n + M p ), g N Nn = -j[g n n-K° + g PP 7T ), Ajv = + V)- 

Furthermore we have introduced the parameter 

< dd > 



(17) 



(18) 



7 



< uu > 



(19) 



to denote the isospin breaking in the quark condensates and set 



1 



< Q1 >= ~[<uu> + < dd>] 



(20) 



From eq.flTB]) we also see that we need to know the value of SX^ 2 , i.e. the isospin breaking 
in the overlaps between the nucleon states and the corresponding interpolating currents. To 
obtain 5Xn 2 , we follow refs. [ffj|9"l,|T0| and use the sum rules for the nucleon two point functions 

(21) 



J d 4 xe ikx (0\T VN (x)rf N (0)\0) — ^Ilf(k 2 ) + 11% \k 2 ) 



which have been considered in the case of isospin breaking in refs. ||15| pTTH . We will take the 
chiral odd sum rule for the amplitudes Ili(A; 2 ) which is known to work better than the chiral 
even ones for n 2 (A; 2 ) fl8] . Including again condensates up to order 4 we have (c.f.eqs.(8) 
and (11) in ref. |7J): 

l/f 2 



, , 4 x 2 V 

(27T — 7- = e^ 1 
4 

v ' 4 



M 6 M 2 g 2 <G 2 > , n , 2 M : 

— + + ( 27r ) — m d <dd> 

M 6 M 2 g 2 <G 2 > , x2 M 2 

— + ^2 + (2n) 2 —m u < uu > 



(22a) 
(22b) 



It should be noted that in this sum rule the gluon condensate g 2 < G 2 > enters in the same 
order as the quark condensate (i.e. order 4) and therefore is taken into account, whereas in 



the sum rule for the three point function (eqs. (|13a|) -(|l~4"|)) it enters two orders higher than 
the quark condensate and was therefore omitted. 

We take the difference between eq. fl22a| ) and eq. (|22 b| ) for p and n and divide by either 
one of them, giving 



5\n^ 
A 



-2) 



'6M; 



N 



M, 



N 



N 



M N J 

2/™ 2 f 2\ 



M 2 



M 2 



M 6 + \g 2 c < G 2 > M 2 



m d -m u 
2 h 7 

m d + m u 



5M N M N 2 
' M N M 2 



(23) 



Putting eq.(|2"3"|) into eq.(|TE|) we obtain the final sum rule 



QNNtt , 



2 4 e 
-3 7 + 37! + 



+ (27^™//* ) 



'5M 



N 



M 



N 



M 2 



2— + 7 - 2 

m d + m u 



SM N M N 2 

Mm M 2 



(24) 



quark mass ratios eivine a value of md . mu 



M 6 + i^ c 2 < G 2 > M 2 
where we have used eq.(|]). 

For the isospin breaking in the quark masses we use the most recent analysis of current 

0.29 ±0.05 

As we can see, one of the crucial ingredients in eq.(^) is the numerical value for the pa- 
rameter 7. Various analyses concerning this quantity have been performed using different 
methods: QCD sum rules for scalar and pseudoscalar mesons [^|.|20|-[22]1 , QCD sum rule anal- 
yses of the the baryon mass splittings [^3| and the D and D* isospin mass differences [24 



as 



well as effective models for QCD incorporating the dynamical breaking of chiral symmetry 
1~6| , |25[| . The range for 7 resulting from these analyses is rather large: 0.002 < —7 < 0.010. 



This range is also consistent with the result obtained from 1-loop chiral perturbation theory 
assuming reasonable values for the strange quark condensate (ss) [Q. 

For the 7r — 77 mixing angle 9 we take the value obtained in lowest order chiral perturbation 
theory ]T4[ : 



1 



-m d - m u 



4 m s — m 



(25) 



Using the numerical values for the quark mass ratios from ref. [[B| we find 9 = (10±0.8)*10 3 . 
Next to leading order corrections are typically of the order 30%, e.g. the decay constants f n 



S 



and f v differ by about 30% if loops are included. It seems therefore appropriate to assign an 
error of 30% to the contribution coming from ir — r\ mixing, i.e. to the term |-^= in eq. ([24]) . 
We have already mentioned that in the treatment of the n-r] mixing we have ignored the 
mixing between r\ and rf as well as n and rj '. The treatment of the rj in the current approach 
is difficult due to the anomaly in the SU(3) singlet pseudoscalar current. The value of the 
n-i] mixing angle 9 increases by about 30%, if the rf is included p6 . 



In order to obtain SMn, we correct the experimental value for the proton and neu- 
tron mass difference by electromagnetic effects, rendering an interval of 1.6MeV < 5Mjy < 



2.4MeV P7U33J|. For g 2 < G 2 > we take the standard value of 0.474GeV 4 , noting that its 



numerical contribution to eq.([24|) is rather small. 

The dashed curve of Fig.|2] shows (~ ~^ J obtained from eq.(|24]) in the Borel window 
0.7GeV 2 < M 2 < 1.5GeV 2 using typical values for the parameters. 

Up to now we have saturated the phenomenological side of the sum rules only with 
the N ground state and have omitted transitions between N and excited N* states as well 
as contributions from the pure continuum. As has been shown e.g. in refs. p8|-ptH] in a 
single variable dispersion sum rule, the transitions N — > N* gives rise to a single pole term 
~ p 2_\ lN i m addition to the double pole term of eq.(^). This single pole term will not be 



suppressed in the Borel sum rules (p4|) . It is easy to see that the inclusion of this contribution 
would add a term to the l.h.s. of eq.([24|) which is of the same general form multiplied by 
an additional power of M , i.e. it can be written as C [j-p) e ~m 7 . The constant C can 
be treated as effective parameter which is optimized in order to obtain the best fit to the 
Borel curve. In the isospin conserving case |7]|||1(| it seems to be justified to neglect this 
contribution due to the fact that the sum rule for gNNn saturated only with the ground 
state is practically independent on the Borel mass M 2 . This indicates that the parameter 
C is compatible with zero. Furthermore the on shell value for g^Nn is reproduced rather 
well in this approach. A recent QCD sum rule analysis for gNNn using two point functions 
PD| also finds that this transition is very small. However, in our case we are looking at 
isospin violation, and there could be a small difference of the parameter C for the proton 
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and neutron contributing to the sum rule (|24]) in the same order of magnitude as ^ . 
It is not difficult to take the excited states into account. The l.h.s. of eq.(pl[) becomes 
(~ 9nn ) +^^ 2 > where the unknown parameter A is optimized in the Borel analysis, which 
means, effectively, by fitting a straight line to the dashed curve of Fig.|^. Doing so results in 
the full line curve of Fig.|] as the final Borel curve for ) ; which is very stable in the 

window under consideration. 

It should be noted that the M 2 dependence of the Borel curve is practically unaffected by 
the large uncertainty in the input parameter 7 and only depends on the ratio , because 

the numerical contribution of 7 as well as 5M^ to the M 2 dependent term in eq.fl24|) is very 
small. The term — 17 + only affects the intersection with the y-axis but not the M 2 
dependence. 

Finally let us look at the effect of a pure continuum starting at a threshold s, which would 
result in multiplying the r.h.s of the eqs.( |13ap -(|HD with the function Ei(x) = 1 — (1 + x)e~ x 
with x = jfj. If one assumes that the continuum thresholds for proton s p and neutron s n 
are equal, there is no effect to the isospin breaking sum rule ([24]). Allowing for a difference of 

= 0.2% (compatible with ^ M -)- l with s n > s p and using a typical value of s = 2.25GeV 2 
would give a contribution of ~ 0.17% to f~ ^~^ ) at M 2 = lGeV 2 . This is noticeably 
smaller than other errors inherent in the sum rule method. 

In order to obtain an estimated error for ( ^— ) we calculate the minimum and the 

V gNNw J 

maximum values obtained from eq. (|24]) after fitting the constant A and using the extreme 
values for the input parameters —7, , 5Mn as well as the the contribution from it — t] 

mixing, as discussed above. This gives an interval of 

17*10" 3 <( ^-)< 30 * 10- 3 . (26) 

Furthermore, from various other isospin violating sum rule analyses (e.g. ref. JT7]]) we know 
that the next higher order condensate (qG ■ aq), which has been omitted here due to the 
reasons mentioned above, may account for about 25% of the leading term. This means 
that we can expect an additional uncertainty of this magnitude. This leaves us with a final 
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interval of 

12 * 1(T 3 < ( — ] < 37 * 1(T 3 . (27) 

The contribution coming from tt-tj mixing, the term in eq.(|24]) amounts to about 8 ± 
2.5 * 10 -3 . As stated above this value would be about 30% larger if 77 — 77' mixing was 
included. The large uncertainty in the input parameter 7 and the lack of phenomenological 
data do not call for a more detailed investigation at the present stage. 

Finally let us compare our result with those of previous studies, which analyze the isospin 
splitting of the pion nucleon couplings arising from the strong interaction, i.e., essentially 
the quark mass difference ma — m u . It should be noted that direct experimental values are 



not available. The Nijmegen phase shift analysis for iViV and iViV scattering data |32,33 



which is consistent with data from nN scattering ||34J| , but includes electromagnetic effects, 



finds ( = 0.002, but with an error of 0.008; thus, there is no evidence for a difference 
and they also find no evidence for a difference between g pnn + and gNN-n- within the statistical 
errors of their analysis. 

From table | we see that our range for f~ ^~ ) i n ( P^D is compatible with the values 
obtained by other authors, both in sign and order of magnitude: 
(1) The quark gluon model of Henley and Zhang [ 35| ; 



(2) the quark pion model of Mitra and Ross |36|,|5|]. This has recently been used by 



Piekarewicz QXJ| , who obtained a violation of the "triangle identity" consistent with the n N 
data analysis of ref. [j3~7f ; 



(3) the use of the quark mass difference — m u and tx — r\ mixing [38[|; 

(4) the chiral bag model [p9 |, which also has our relation (eq.(|15|)) for the charged coupling 
to be valid. 



(5) On the other hand, the use of the cloudy bag model |4(J leads to J ~ 0.006, with 

the opposite sign to our result. 



It should be noted that there are electromagnetic corrections, whose direction are un- 
known. The charge difference we obtain due to the strong interaction would, by itself, lead 
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to a difference on the scattering lengths |o nn | — \a pp \ ~ —0.5 ± 0.2fm, smaller than, but in 
the opposite direction to the observed difference ||. Of course, there are other effects which 
play a role, e.g. p — uj mixing. 

We are aware that using the three-point function is a priori less suitable than the two- 
point function for calculating the pion nucleon coupling on shell, because it works at spacelike 
q 2 ~ — lGeV 2 and needs the detour of comparing the \ pole residua |]7|j9|Jl0l| . As we have 
already mentioned earlier this means that one has to assume that n pole dominance can 
still be applied in this region and higher pseudoscalar resonances are neglected, or, in other 
words we use the PCAC interpolating pseudoscalar field at those values of q 2 . A quantitative 
analysis of the contribution of these higher resonances would require some knowledge about 
their coupling to the nucleon. There are various indications and consistency checks that the 
concept proposed in refs. |]7|j9|JT0]| is a reasonable one: First of all in the isospin violating case 
one is likely to be less sensitive to the higher resonance contributions due to cancellation 
which presumably occur if summing over the higher pionic excitations. Furthermore in 
the isospin conserving case JT|,|9|,[iTi|] the experimental value of gNNw is reproduced rather 



well. The Borel stability in both the isospin conserving J7^[5|,rO| and the isospin violating 
case (this work) is a further consistency check, although this of course only a necessary 
but not sufficient condition. Finally there is the analysis of the ttNN formfactor |y| using 
this approach. The q 2 dependence of gNNnil 2 ), which contains effectively the higher pionic 
resonances in the spectral function ([!]), is consistent with various other approaches using the 
same interpolating current but working at lower q 2 . This indicates that the interpolation 
between low and high q 2 region is done reasonably well. Despite all these arguments the 
importance of the higher pseudoscalar resonances remains a matter to be settled and needs 
further quantitative investigation pE2f . 

On the other hand the use of a two-point function [|7|, p!0| , ^T| , |30| has other, and we believe 



worse problems in the consideration of isospin violation. The single nucleon pole sum rule, as 
it has been used in refs. []T|, |TU| , pET| , suffers principally from the problem that the contribution 



from the transition iV — > N* enters exactly in the same form as gNNn itself, namely as 
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single pole. After Borel transform one obtains a term gNNn + A instead of gNNn + AM 2 
as in case of the double pole sum rule. Hence within the single pole sum rule itself there 
is a priori no way to separate the N — > TV* contribution A from gNNn- In refs. |T|, |iO|j4T 



the N — > N* transition has been ignored. In the isospin conserving case this is a posteriori 
justified because the numerical value of this term turns out to be small, as it has been 
discussed above. However we do not know if this is true in the isospin violating case. The 



double nucleon pole sum rule |30j would avoid this problem and moreover is able to give 
a value for the iV — > N* contribution. Unfortunately as one can see from the analysis in 
ref. 0] already in the isospin conserving case this sum rule seems to be rather sensitive 
to the condensate input, which is the quark condensate < qq > and especially the higher 
order mixed condensate (OlqGj^qln) , which has to be included to give a reasonable value for 
gNN-w- In case of the three point function the condensate input parameters are much better 
under control. For these reasons we prefer to work with the three point function sum rule. 

To summarize we have calculated the splitting between the pion nucleon coupling con- 
stants g ppn o, g nri7r o and g vn - n + due to isospin breaking in the strong interaction by using the 
QCD sum rules for the corresponding pion nucleon three point functions. We have taken 
OPE diagrams up to order 4 into account. Our result for the splitting in the neutral cou- 



j-2 ^ 9 PP tt° 9 nn7T ^ q 7^ m-2 

gNNn 

average of the two neutral ones. 



plings is 1.2 * 10 2 < pp7r gNN nn7T < 3.7* 10 2 . The charged coupling g pnn + is found to be the 
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FIGURES 

FIG. 1. Diagrams in the OPE. 

FIG. 2. Dependence of {g pp -K — gnmr ) / 'ONNw on the square of the Borel mass M 2 . As example 
we have used the parameters 7 = -0.01, 5M N = 2 MeV and = 0.28. The dashed curve is 

obtained by omitting the transitions N — > iV*. In the full curve these contributions are included. 
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TABLES 

TABLE I. Comparison of (g P pn ~ 9nnn ) I '9NNn obtained in different approaches including 
isospin violation effects from strong interaction. In order to compare the numerical values of 
refs. fl|,|35| with the other results we have used SMn/Mn = 0.002. 

)/9NNtt 

this work ss 0.012 . . . 0.037 



Ref. |35| 0.010... 0.014 

Ref. |H w 0.006 

Ref. || 0.005 ± 0.0018 
Ref. M » 0.0067 
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